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Abstract: The Dirac-Born-Infeld (DBI) action has been widely studied as an interesting 
example of a model of k-inflation in which the sound speed of the cosmological perturbations 
differs from unity. In this article we consider a scalar-tensor theory in which the matter 
component is a field with a DBI action. Transforming to the Einstein frame, we explore 
the effect of the resulting coupling on the background dynamics of the fields and the first- 
order perturbations. We find that the coupling forces the scalar field into the minimum 
of its effective potential. While the additional scalar field contributes significantly to the 
energy density during inflation, the dynamics are determined by the DBI field, which has 
the interesting effect of increasing the number of efolds of inflation and decreasing the 
boost factor of the DBI field. Focusing on this case, we show, with the benefit of numerical 
examples, that the power spectrum of the primordial perturbations is determined by the 
behaviour of the perturbations of the modified DBI field. 
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1. Introduction 

The popularity of string/M theory and its interpretation as a truly fundamental model 
of physical phenomena has motivated the search for a mechanism to explain the process 
of inflation. In canonical models of inflation, a slowly rolling scalar field, with its energy 
density dominated by its potential, is introduced as a phenomenological tool to account for 
the dramatic expansion of the universe to solve the flatness and horizon problems of the Big 
Bang theory and also as a source of the primordial fluctuations that give rise to structure 
formation and the observed anisotropies in the CMB radiation. While there is no shortage 
of scalar fields in stringy models, problems arise when one looks for a suitable candidate 
to play the role of inflaton in the canonical models |]l|, |2| . The abundance of light scalars 
(moduli) not only complicates the dynamics of inflation, but also means that models must 
be tuned to prevent unwanted light moduli affecting the post-inflationary universe. As well 
as this, although heavy fields are effectively frozen during inflation they may contribute to 
the potential that determines the evolution of the dynamical fields, further complicating 
the 'eta problem' in inflationary model building, in which one encounters large corrections 
to the flat potential required in slow-roll inflation 1^, ^, |5|, |6|, 0. 

An intriguing scenario is that of Dirac-Born-Infeld (DBI) inflation |^, |9| in which the 
role of inflaton is played by the radial coordinate of a D3 brane moving in a warped region 
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(throat) of a compactified space. Crucially, a speed limit is imposed upon of the motion 
of the brane (which appears as a point-like object in the compact dimensions), affected by 
both its speed and the warp factor of the throat. The speed limit introduces a parameter 7, 
which plays an analogous role to the Lorentz factor in special relativity, that grows without 
bound as the speed limit is approached. This has the effect that in the 'relativistic' regime 
(as we shall call it), the importance of the friction and potential terms in the scalar field 
equation of motion is lessened, so that a DBI field can drive inflation with a relatively steep 
potential. 

From a cosmological point of view, DBI inflation is important as it is a concrete 
example of k-inflation i.e. inflation driven by a scalar field with a non-canonical kinetic 
term 0, |ll|, ||, ||, . A distinguishing feature of these models is that the 

field perturbations travel at a speed less than that of light, with the result that the Fourier 
modes freeze in at a sound horizon that is different from the curvature radius and which 
changes dynamically with the background evolution This is an exciting possibility, 

as the amplitude of the primordial perturbations would not be (almost) independent of 
wavenumber as in the standard inflationary models and the effects of the evolution of the 
universe during the inflationary regime could be encoded in mode correlations. In DBI 
inflation, the sound speed is given by 7"^, and decreases to small values as the speed limit 
is saturated. Observational limits on the level of non-Gaussianities in the CMB radiation 
can thus be combined with traditional constraints arising from the power spectrum of 
primordial perturbations to distinguish between DBI and canonical models of inflation 
(see ^ for recent reviews). 

It is to be expected that the effective action governing the inflationary phenomenon will 
contain other features characteristic of fundamental theories such as non-minimal coupling 
and multiple fields. In particular, the effect of the presence of many scalar fields on the 
inflationary dynamics has been the subject of intense study over the past decade, with 
many works focusing on the behaviour of entropic or isocurvature perturbations in these 
models. These can source the adiabatic fluctuations in the system so that the curvature 
perturbation, which in single field models is conserved outside the horizon, acquires a 
dependence on the inflationary dynamics. Interestingly, the presence of a coupling of the 
type A{<j)) = exp(/?(/)) between two scalar fields (arising naturally in scalar-tensor theories 
considered in the Einstein frame) has been found to further influence the transfer between 
entropic modes and adiabatic modes |21, 22|; however, the analogous case for multifield 
k-inflation models has been neglected in the literature. Multiple fields arise naturally 
in DBI scenarios, where the extra degrees of freedom can be identified with the non- 
radial coordinates of the D3 brane in the compactified space, and several studies (both 
analytical and numerical) have been devoted to understanding the evolution of adiabatic 
and isocurvature perturbations and the predicted level of non-Gaussianity in these models 
|23, 24, 25, 26, 27, 29, In contrast, the question of the behaviour of the perturbations 



in models involving many fields, each with a different intrinsic sound speed, has received 
relatively little attention. One such model was investigated in [31|, in which inflation is 
driven by two standard DBI fields, each corresponding to a brane with its own sound 
speed dependent on its evolution. The presence of two sound speeds means that the 
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perturbations freeze in at the sound horizon corresponding to the larger value of Cs while 
the non-Gaussianities are at a level corresponding to the square of the larger value of 7. 

In order to address these issues, we propose a model in which a scalar field with 
a DBI action evolves in the presence of a canonical scalar field. We embed the DBI 
inflationary scenario into a scalar-tensor theory, in which the additional field could describe 
the degree of freedom associated with additional moduli fields in the higher-dimensional 
theory. We will not construct such a model from first principles, but rather treat it as a 
phenomenological model. Therefore, we consider the generalised case in which the fields 
are coupled, motivated by the possibility of the additional scalar field having a non-minimal 
coupling to gravity in the effective action. Notwithstanding the importance of non-minimal 
coupling from a fundamental theory standpoint, couplings of this type have often been 
neglected in the DBI inflation literature, partly due to the success of the standard DBI 
scenario as an consistent phenomenological model in their absence. However, we should 
point out that there have been some interesting studies such as |p, where non-minimally 
coupled DBI models were considered in order to circumvent the eta problem in inflationary 
model building. Some issues pertaining to the realisation of non-minimal coupling in DBI 
models of inflation from the string theory perspective were raised in and also . The 
authors of performed a detailed analysis of a single field DBI scenario with non-minimal 
coupling of the Brans-Dicke type. Their analytical work shows that both in the small 7 
and large 7 limits the effect of the coupling is negligible but in the intermediate regime, 
the degeneracy of observables can be broken: for a given tensor-to-scalar ratio and medium 
/nl, non-minimally coupled DBI models have a redder scalar spectral index. The effect 
of coupling in multifield DBI models was also considered in |34, 35 1, focusing on the effect 
of particle production due the interaction between the 'inflaton' brane and trapped branes 
in the warped throat. Although, like our model, this scenario involves a DBI action non- 
minimally coupled to a second field, the motivation and dynamics are very different: for 
example, in our model the scalar field as well as the DBI field can contribute significantly 
to the energy density and the background dynamics depend strongly on the evolution of its 
vacuum expectation value (VEV). Recently, coupled DBI infiation has also been considered 
in the context of a warm inflationary model | |36| . 

This paper is organised as follows. In Sec. ^ we introduce the model in the Jordan and 
Einstein frames, and derive the equations of motion for the inflationary Universe. Sec. ^ 
is devoted to an analysis of the background dynamics, both in the non-minimally coupled 
and coupled cases and in Sec. § the equations determining the evolution of the first-order 
perturbations are derived and solved numerically. We conclude in Sec. ^ with a short 
discussion. 



2. The Model 

Non-minimally coupled scalar fields arise in a number of scenarios in high-energy physics, 
such as low-energy effective actions from higher dimensional theories, quantum field theory 
in curved space and f{R) models of gravity. We consider a scalar-tensor theory with 
canonical field if coupled to the Ricci scalar combined with a second field x with a DBI 
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action. In the Jordan frame, this is described by the action 

S = [ d^x^g \f{^)^ - \{V^f - U{^) + r\x) [1 - 7^'] - y{x) 



(2.1) 



where the function F{'^) determines the coupHng of the canonical scalar field to gravity, 
f{x) is the warp factor of the DBI field, U{ip) and V{x) ^ire potentials for the canonical 
and DBI fields respectively and 7 is the boost factor of the DBI field, given by 



1 



7 



Vi + /(x)(Vx)2 

Varying ( p.l| ) with respect to the metric gives the modified Einstein equations. 



(2.2) 



FG^^ = T°J^i + [d^^d^^ - \g^^{V^f\ - gf,^U + V^^V^F - g/^^aF, (2.3) 

where □ = g^'^V^V^ is the d'Alembertian operator and T^^^ the energy momentum tensor 
of the DBI field. The equations of motion for the fields are 



2-cuaip = F'T - w'iViff - mF' + 2U'F , 



(2.4) 
(2.5) 



where zu = F + |F'^. In the Jordan Frame, the equations governing the dynamics are 
evidently quite complicated. In order to simplify the system and gain a different perspective 
we can perform a conformal transformation to the Einstein frame in which the Hilbert term 
in the action takes the standard form. We define 



(2.6) 



where the tilde denotes quantities evaluated in the new frame in which the field is 
minimally coupled to gravity. Using the definitions. 



d(p 
dip 



F, 
F 



+ 



U{if) = U{^)F-\ 
^(^) = F-5, 
1 



7 



the action in the Einstein frame can be written as 



S 



d^x 



where, 



(2.7) 

(2.8) 
(2.9) 

(2.10) 



(2.11) 



Sdbi [aV] = / d'xy^A^ {rHx) (1 - 7"') - VU)} . (2.12) 
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The scalar field can thus be redefined to retain its canonical kinetic term in this frame; the 
DBI field is not redefined but its behaviour will be affected by the coupling terms A((p) 
present both as an overall multiplier in the action but also under the square root in the 
definition of the boost factor 7. Throughout the rest of this paper we shall work in the 
Einstein frame and so will neglect the tildes for convenience of notation. 

We study the fiat FRW metric in this frame, describing an expanding universe, 

ds"^ = -df + a^{t)5ijdx'dx^ , (2.13) 

with scale factor a{t). The background equation of motion for the homogeneous field x iS) 

X + m^-^x + - 37"' + 27^') + A^r^V^ = -/3X9i(37"' - 1), (2.14) 

where dots indicate derivatives with respect to the time coordinate in the Einstein frame, 
= df /dx, H = a/a and /3 = dlnA/dip. Throughout this paper we will assume that /3 is 
constant, so the coupling A{ip) is given by 

A{ip) = eMM- (2.15) 

Eqn. ( 2.14| ) reduces to the standard DBI equation of motion when /3 = 0, although it is 
important to note that in the coupled case, the boost factor 7 acquires a factor of A under 
the square root that affects the dynamics of the field, 

^ (2.16) 



The (f equation takes the simple form, 

¥3 + + [/^ = /^Tdbi, (2.17) 

in terms of the trace of the stress-energy tensor of the DBI field Tdbi = — Pdbi + 3pDBi • 
The energy density and pressure for x c^re 

PDBi = A^[f-\j-l) + V] PBBi = A^[f-\l-j-^)-V]. (2.18) 

so Tdbi can be written in terms of 7, A and V as 

Tdbi = [/-'(4 - 87"' - 7) - 4F] . (2.19) 

The Priedmann equations are 

3H^ = ^if^ + U + PBBI, (2.20) 

-2H = ip'^ + -fA^X^. (2.21) 

From these we can define the slow-roll parameter, e, 

_ H 3 ip^ + -iA^x^ 



2pDBI + ^V^2 + ^' 

in the standard way, so that d/a = H^{1 — e). 



(2.22) 
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Figure 1: DBI inflation, plotted against iV = Ina. The upper left panel shows the log of the field 
value of x; the upper right panel shows the increase in the boost factor 7. The lower left panel 
shows the logarithm of the slow-roll parameter e, and the bottom right shows the equation of state 
parameter wdbi = vl P-, which remains close to —1 until x starts to oscillate about its minimum. 
A reasonable number of cfolds of inflation are generated in this case, but the boost factor gets very 
large. Parameter Values: A = 10^^, m = 5 x 10^^, /i = 0.1. 



3. Background Dynamics 
3.1 DBI inflation 

Before discussing the two-field system, it will be useful to remind tlie reader of some of the 
features of single- field DBI inflation HI . Setting /3 = in ( |2.14 ) the equation of motion for 
the minimally coupled DBI field can be written as 



X + mr'^X + - 37"' + 27"') + 7"'^x = 0. (3.1) 



As the X field rolls down its potential, 7 gradually increases from 1, driven by the potential. 
As Ixl increases, so does the boost factor, with the result that terms in I/7 are suppressed. 
The potential and friction terms then drop out and the equation reduces to. 



X + i||«0 ^ X - ^ (x' + Z^'x) ~ 0, (3.2) 



where in the second equation we have used the standard 'mass gap' form for the warp 
factor, given by 

This is an approximation to the Klebanov-Strassler throat (see [^ ] and references therein) 
used to describe the warped conifold in which the DBI field moves that reproduces phe- 
nomenologically the form of /(x) in both the IR and UV limits. In the limit of an AdS 
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throat {x^ 3> //^) the solution is x ~ \/A/t; towards the tip of the throat where 
X decreases faster as x ~ exp(— /iy^2/At). 
The slow-roll parameter e is given by, 

2/-i(7-l) + 1^ 27 + /F ^ ^ 

Note that although 7 gets large, the large value of / required to satisfy observational 
constraints Q means that the potential term is the dominant contribution to the energy 
density and the pressure, implying that e is less than 1 when x is small and the equation 
of state u;dbi = Pdbi/pdbi 

™ = f-\i-i) + v " -^TTF ^^-'^ 

becomes very close to —1. This is plotted in Fig. ||. 

Single field DBI inflation is normally treated analytically using the Hamilton-Jacobi 
formalism, in which the field value x assumes the role of the time-coordinate. Differentiat- 
ing H = H{x) with respect to cosmic time and using the second Friedmann equation, one 
finds, 

*^-^w <^-^' 

Therefore, the definition of 7 can be rewritten as. 



lix) = Vr+47(xWxF, (3.7) 

and the potential, 

v{x) = mix? - f{xr\i{x) - 1). (3.8) 

The slow-roll parameter ^ can then be written 

2 fH'V fs 1 



where the second equality holds for the AdS throat with potential w?x^ when e ^ 1. An 
unusual feature of DBI inflation is that an increased mass term decreases e, unlike in the 
standard slow-roll case. The advantage of the Hamilton-Jacobi approach is that, given a 
potential ^(x), eqns. ( ^^ ) (using the definition of 7(x) in ( |3.7D ) and ( |3.6| ) can be solved 
sequentially. One can also choose an anzatz such as oc x in order to calculate V . 

In much of the DBI literature w?x^ potentials similar to those in large field inflation 
are used, following the analysis in P]. This showed not only could a mass term be generated 



^One can generally construct a 'tower' of flow (slow-roll) parameters (the first two being e and r\) defined 
in terms of derivatives of i?(x) with respect to x §]• In addition, it has been shown in that one can 
additionally construct a second 'tower' in terms of derivatives of 7(x) with respect to x- The first flow 
parameter in this 'tower', s, enters the expression for the spectral index of curvature perturbations as a 
measure of the variation of the speed of sound i.e. Us — \ = — 4e -\-2r] — 2s. 
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in the DBI scenario without the setup being ruined by back-reaction of the probe brane, 
but it could drive the power-law inflation discussed above. Although it has since been 



found 1 40, 41, 42| that the class of potentials that lead to successful DBI inflation is much 
wider, in this work we focus on modifications to the standard DBI scenario and will restrict 
ourselves to the potential 

V = \m\\ (3.10) 

In the cutoff throat, the 7 increases to a maximum then decreases slightly as the friction 
term starts to become more important. Generally, 7 remains large even after x has passed 
through the origin, whereupon the other terms in the equation start to combine to stop the 
field, giving 7 = 1. After this point, the field exhibits oscillations about its minimum. The 
presence of one or more D3 branes in the tip region of many throat constructions, means 
that DBI inflation is often taken to end with the probe brane annihilating at the point 
Xend) rather than with a violation of the condition for acceleration e < 1. The mechanism 
by which the energy released in this scenario is transferred to the standard model brane is 
discussed in ^3|. If x can reach the bottom of its potential, it can oscillate, which allows 
for the possibility of energy transfer in a period of (p)reheating p^ , however, this is 
beyond the scope of the present paper. 

As discussed in the introduction, DBI inflation is a particular case of k-inflation in 
which the squared sound speed of scalar perturbations is equal to = 7"^. The non- 
Gaussianity of the primordial perturbations is determined by the bispectrum (the three- 
point correlation function of the scalar curvature perturbation C), which in general is a 
non-trivial function of the momentum configuration. The bispectrum can be decomposed 
into several parts (shapes), each with a magnitude parameterised by an estimator (or non- 
linearity parameter) /at^ defined by an ansatz of the form [^] 

c = a - \fNLCh (3.11) 
5 

where is the linear part of the scalar curvature perturbation. In DBI inflation, the 
leading contribution is estimated by^ 

35/1 
108 



fNL = — [-,-l]. (3.12) 



Therefore, 

fNL ~ 0.3272. (3.13) 

Conservative estimates of the observational value^ of the non-linearity parameter have 
I/a^-lI ^ 0(100), leading to an upper bound on the value of the boost factor 7 < 20, 



^Eqn. ( 3.12 ) is calculated on an equilateral triangle using momenta of equal amplitude. This is relevant 
because the scalar perturbations are frozen on large scales and oscillate (averaging out) on small scales. 
Thus, the dominant contribution comes from modes with similar wavelengths that exit the (sound) horizon 
at the same time. 

^We note in passing that the measured value of /jvl can from receive relevant corrections from the 
post-inflationary evolution of the perturbations and differ from the expressions given above. [^. The 
curvature perturbation, however, is conserved on large scales so one can connect the observable signal to 
the primordial perturbations. 
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Figure 2: Minimal coupling with a quadratic potential U{(p) = Uq(p'^ with Uq = 10~^^ (black, 
dotted), 10~^^ (blue, dot-dashed), 2 x 10~^^ (red, dashed) and 5 x 10"-^° (cyan, solid). The upper 
panel shows the boost factor 7, the middle panel show the DBI density parameter r^DBi = /0DBi/3i?^ 
and the lower panel the slow-roll parameter e. Plots are against efold number N — ln(max(a)/a) 
so = corresponds to the end of inflation. Other parameter values: A = 2 x fO^^, m = 5 x fO~^ 
and II = 0.2. 



assumed to apply during the first 10 of the last 60 efolds of inflation which directly affect 
the CMB radiation. Since the boost factor increases with the number of efolds of inflation 

sufficiently large predict a high level of non-Gaussianity 
I M and also |, pi, |3l). 



A'^max, many DBI models with Nj^ 
exceeding current observational bounds (cf. 



3.2 Minimal Coupling (3 = 

Let us now include the field. A priori, there appear to be five possible ways in which 
the system can behave, viz. x dominated inflation, (p dominated inflation, a period of 
double inflation with a changeover from x to </j (or vice-versa) or a setup in which x a-nd 
ip contribute roughly equally to the energy density. 

In the minimally coupled case, the fields x ^ only interact gravitationally and 
since they evolve on different time scales in this case, it is difficult to get a situation in 
which both contribute equally to the energy density for a considerable length of time. Even 
if one fine-tunes the initial conditions so the fields have similar energy densities, there are 
problems, x becomes small quickly unless its potential is large, but if we require 
(where ilj = pi/3H^ is the density parameter for the i*^ field) the potential for ip is steep 
and the field quickly rolls to its minimum. The opposite is true if y ~ [/ is smaller: since 
the DBI contribution to the slow-roll parameter depends explicitly on the DBI field mass, 
it decays quickly, leaving ip to slowly roll to its minimum. 

If the total energy density is dominated by the DBI potential, the evolution of x 
proceeds in the usual fashion: 7 increases steadily to large values, until X A*! whereupon 
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2 4 6 8 10 

Figure 3: The evolution of thie canonical and DBI fields in the ((/?, x) plane for the plots shown in 
fig. i 

the field exhibits oscillations. If the difference in energy densities between the two fields in 
not great, the field rolls in the background. If, however, 2>H^ ~ ^(x)) Hubble damping of 
the ip field is disproportionally large, so the field evolves slowly until the DBI field decays. 
When this occurs, the ip field rolls to its minimum as in canonical slow-roll inflation. An 
example of this behaviour is the black (dotted) curve in figs. I and |. 

The opposite occurs when the ip field is initially dominant (i.e. , ~ U{ip)): one 
finds a period of standard slow-roll expansion followed by DBI infiation (cf. solid cyan curve 
in figs. |2| and ^). At first O,^ > and 7 ~ 1 so the DBI field is dynamically unimportant, 
simply frozen or evolving slowly to smaller values. The timescale for ip infiation is rapid 
compared to that of DBI infiation so 7 does not increase much during this stage if x is 
rolling. When the <p field decays, the slow-roll parameter grows, only to decrease again as 
the boost factor increases. In both cases, if one field is dominant, there are two distinct 
periods of infiation. 

Perturbations in two-field infiationary models can often be decomposed into modes 
tangential (adiabatic) and orthogonal (entropic) to the background trajectory in field space 



-10- 



|5C, 51, |22|. The trajectories for the two aforementioned examples exhibit a sharp turn 
in field space (fig. |3|) as the dominant field reaches its minimum. In canonical scalar 
field models, behaviour of this type is of particular interest as it can lead to a transfer 
between entropic (isocurvature) and adiabatic modes, so that the curvature perturbation 
is not constant on super-Hubble scales. If one fine-tunes the mass scales so that both 
fields contribute roughly equally to the total energy density for a considerable number 
of efolds, the dynamics of the DBI field can be significantly affected by the additional 
contribution to the Hubble damping, leading to smaller values of 7 (relative to the DBI 
dominant case). As both fields are evolving, the interplay between the field perturbations 
suggests that isocurvature fluctuations could be important. However, minimally coupled 
case is extremely sensitive to small changes in the field masses (as can be seen in fig. ^ 
and the initial conditions of both fields so will be analysed in more detail in future work. 
In the rest of this paper we will focus on the non-minimally coupled case, in which the ip 
field moves in an effective potential dependent on the DBI field. 

3.3 Non-minimal Coupling 

If we consider the non-minimally coupled model interactions between the fields can become 
more complex as the extra mass term — /STdbi in ( ^-17 ) starts to play an important role. 



The combined effect of this and the bare potential means that the effective potential C/cff 
U — ^Tdbi can have a minimum value at which the field value ip = (pmm satisfies. 



dU 

dtp 



_ P^^PVrmnrpb^^^ ^ (3.14) 

'/5=</'min 



where T^^j = /~"^(4 — 7 — 37^"*^) — AV ~ —4V is the trace of the 'bare' DBI stress-energy 
tensor and the ip dependence of A{(p) has been written explicitly using ( ^.151) . 

The existence of a minimum is dependent on the form of the bare potential: the 
exponential form of the coupling means that it is necessary to have a potential that slopes 
to the right^ i.e dU/dp < 0. 

3.3.1 Offset quadratic potential 

In this paper, we will consider two choices of potential, dependent on two free parameters 
Uq and r]. First, the offset quadratic potential 

U{ip) = Uo{p-v)\ (3.15) 

for which the true minimum is p = ij. (The advantages of the offset potential over the 
standard quadratic r] = case will become apparent.) In this case dU/dp = 2Uq{p — r]); 
writing v = —4:/3pmin, the condition for the minimum is 

,__2^^7jBie^ 8/3^y(x)e^^^ 



^Alternatively, if 13 were defined to be negative, for a minimum to exist one would need a potential with 
positive slope dU /dip < 0. In this case, one would need negative values of ip to ensure that A > 1, but the 
phenomenology of the model would be unchanged. 
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Figure 4: Non-minimal coupling with yl < 1 for the quadratic potential. Black dot-dash curves 
have (3 — 0, blue dashed curves have /3 — 0.25 and solid red curves have (3 — 0.5. Note that in the 
latter case the effective potential is steep and there is a brief kination stage during which there is 
no inflation (e > 1) as is large. Parameter values: A = 2 x 10^^, m = 5 x 10~^, Uq = 2 x 10~^^, 
^ = 0.2, 7? = 0. 



where the approximation is valid as, generally, the pressure and energy density of the DBI 
field are dominated by its potential. ^Pmin can then be expressed in terms of Lambert's W 
function W{x) (satisfying W{x)e^^^^ = x) as 

The effect of the coupling on the dynamics of the system depends strongly on the sign 
of 99, as this affects whether the coupling A = exp(Pip) is greater or less than one. Even if 
ip is initially positive, the field will be attracted to its minimum, which can take negative 
values if rj is small. As can be seen in fig. ^, while ip is being dragged to smaller values the 
extra friction term on the RHS of ( ^.14 ), which is approximately f3(fX7 ^-cts to further slow 
the rate of decrease of x- When if reaches its minimum, depending on the steepness of the 
effective potential it will either overshoot (which in extreme cases can lead to values of A 
substantially less than one) or track the minimum value, which evolves slowly toward r] at 
a rate determined by the decrease in V{x)- 

In the non-minimally coupled case, we find numerically that for a wide range of param- 
eters if quickly settles in the minimum of its effective potential, even when substantially 
displaced (cf. fig. This considerably simplifies the analysis as one can make the as- 
sumption (p = V9min for most of the time during inflation, so the background dynamics are 
largely independent of the initial conditions for the scalar field (p. To better understand 
the dynamics in this case, consider the condition for the minimum ( p. 14 ) for this potential, 

2Uo{ip - 7?) ~ -4/3^ V. (3.17) 
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Using this, we can rewrite the potentials in terms of W = W{x) with x 



8^2 



U 



I6/32 



so the Hubble parameter becomes 



m"^ ~ (1 + \W)AW. 
One can differentiate ( 3.17| ) with respect to time to get 99 in terms of x 

J_ ( W \ X 

2/3 \1 + W) X 



8^2y^iPv/Uo as 
(3.18) 

(3.19) 

(3.20) 



where we have used the quadratic form of the x potential to write V^/V = 2x- Thus, the 
rate of change of the scalar field ip has an explicit dependence on x- Focusing on the rate of 
change of the Hubble parameter, it can be seen from eqn. (2.21) that for the kinetic term of 
the DBI field to dominate over that of the scalar field, we need the condition A'^jx^ ^ 992 
to be satisfied i.e. 

2 



^ 7X > 



7 > 



7 > 



1 



W 



l + W 

w 



XT 

X^ 
2 



?2 \ l + W 

(2 + W)W^ 
48/32(1 + Ty)2 



1 



A^V 
A^m^ 
IP' 



(3.21) 



where ( p.l9| ) was used in the third step. Notice that since VF/4/3 is equal to the dislocation 
of <p from its minimum, W/2 is fairly small, around 0(1). Since the effective mass of the 
DBI field is mA [cf. ( p. 14 )], which should be (at least) the order of the Hubble scale, we 
find that the ip field is dynamically subdominant (i.e. -H = ^(99^ + A?^x^) « A^^x^) if 
7 ^ C(l)- So, when x is affected by the DBI terms in the action, the slow-roll parameter 
becomes 



1 



37X" 
2V 



The case in which 



{2 + W)A^V ^2(1 + 1^) 
> is thus of particular interest as this gives A> 1, yielding 
'37X'" 



(3.22) 



< 



2V 



where eDBi is the slow-roll parameter for a standard DBI field^ of mass m. 
inversely related to the number of efolds of inflation A^max = ln(aj/aj) by 



Since e is 



N, 



max 



f' Hdt= [ ^ -^\/^dtc^ r"'' Ay 
Jti Jti Jo 



-dx, 



^It should be noted that the factor oi A ^ in the definition of 7 in the Einstein frame means that 7 is 
smaller in the coupled case compared to the standard DBI scenario so even discounting the small factor 
A-^{1 + in (I^H), ^'yx'^/V is generally less than eDBi- 
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Figure 5: The effect of increasing /3 from 0.1 to 1.5 for the offset potential case. The upper-left 
panel shows the trajectory in field space (starting at % = 1.5) with tp tracking the minimum, the 
upper-right panel shows the DBI boost factor 7, the lower- left panel shows the slow- roll parameter 
rj and the lower-right panel shows the density parameter of the DBI field. It can be seen that 
increasing the coupling increases the duration of inflation and decreases the maximum boost factor. 
Parameter values: A = 8 x lO^^^ m = 8 x 10'^, Uq = 10"^°, = 0.15, = 3. 



this results in a increase in the number of efolds of inflation. This is evident in fig. |5|, in 
which the the evolution of the fields is plotted for several values of /3. One can see that 
increasing the coupling from 0.1 to 1.5 decreases the slow-roll parameter by over an order 
of magnitude, causing the the number of efolds of inflation to increase from ~ 30 to ~ 140. 
Also, the maximum value of the DBI boost factor is dramatically reduced by a factor of 
10, which could lead to interesting observational deviations from the standard DBI model 
(discussed in the final section). 

If ^ <C 1, the effective mass of the DBI field is decreased, which has the twofold effect 
of increasing the slow-roll parameter and decreasing the number of efolds of inflation, as 
well as affecting the DBI speed limit via the factor of A"'^ under the square root in eqn. 
( 2.16 ). This can be seen in the cases presented in fig. ^. As ifmin is negative, ^ ^ 1 and 
the number of efolds of inflation is decreased compared to the minimally coupled case. 

One can then ask, what conditions need to be fulfilled for in order for the coupling A 
to be greater than one? Substituting (|3.16| ) into (pmin > gives the condition 



W < 4/3??, 



(3.23) 



which is plotted in the left-panel of fig. ^ for several values of the ratio V/Uo. As one might 
expect, it is easier to satisfy this condition with larger values^ of the bare minimum -q, as 



®Note that the mexican hat potential with minima at positive and negative values of (p is not suitable 
here as, if the field starts at small values, it is dragged in the direction of the negative minimum by the 



term on the RHS of eqn. (2.17) giving A < 1 
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Figure 6: Conditions for a positive minimum and DBI dominance for the offset quadratic (left 
panel) and exponential (right panel) potentials. Left panel: regions above the solid lines satisfy 
the condition for i^rnm > (3.2S) using V/Uq = 0.01, 0.1, 1, 10 (starting from the bottom). Shaded 



regions indicate where the DBI energy is dominant (3.24) with the lightest region corresponding 
to V/Uq = 0.01. Right panel: the equivalent conditions for the exponential potential. In this case 
the region satisfying the condition for DBI dominance ( |3.30D , represented by the shaded region, is 
independent of the ratio V/Uq. 



the minimum of the effective potential is shifted to larger (more positive) values. Using 
( 3.19 ), the condition for DBI dominance, J^dbi > is 



W <2, (3.24) 

which is also plotted in fig. ^. Naively, one might expect the situation wherein the DBI 
plays a dominant role to be characterised by V/Uq ^ 1. However, this means that the 
second term in ( 3.16| ) is dominant unless the offset is large so the net effect is ipmin < 0, 



leading to smaller values of the energy density ~ A^{ip)V. One can see in fig. I that 
the region of the P,r] parameter space that can satisfy both ( p.23 ) and ( 3.24| ) is fairly 



small meaning that for a large number of cases for which ^4 > 1 the DBI energy density 
is subdominant. However, this does not mean that x is dynamically unimportant: the 
opposite is the case. The large contribution of the energy density of the scalar field arises 
as cp is trapped in a false vacuum away from its true minimum, so if cp tracks the minimum, 
its evolution is determined solely by the DBI field. This can be seen in the examples shown 
in fig. ^ In each of the cases presented here r^DBi ^ 1, but the field is never completely 
subdominant. However, despite the significant contribution of il^ to the energy budget, 
the dynamics are controlled by the DBI field and if is trapped in the ^-dependent minimum 
of its effective potential. (This behaviour is similar to hybrid inflation [^] , see also [^ ] for 
a recent discussion of hybrid inflation in scalar-tensor theory.) 
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3.3.2 Exponential potential 

As well as the offset quadratic potential, we also consider an exponential potential 

Ui^) = Uoexp{-r]ip), (3.25) 

in which the parameter rj determines the slope and, under the inffuence of the bare potential, 
ip moves to larger values. The condition for the minimum in this case reads 

U ~ ^A'^V, (3.26) 
V 

so the minimum is given by the logarithmic function 

Note that in this the fields are dominated by their potentials, the energy densities 

of the fields are proportional to each other. Thus, we have a tracker solution for the fields, 
in which the density parameters are constant f^DBi oc $7^ ~ const (as can be seen in fig. 
0). The condition for 9?min > is therefore 

r] > 4./3{V/Uo). (3.28) 

This is plotted in the right-panel of fig. ^, and one can see that the condition is similar to 
that for the offset quadratic potential, where rj was the position of the true minimum. This 
is simply because in both cases increasing the parameter r] steepens the bare potential at 
small values of (p. 

Using ( 3.17 ), the Friedmann equation can be written, 

3^2 = (l + 4/3/r/)^V, (3.29) 

so the condition for DBI dominance is independent of the ratio V/Uq: 

r] > 4/3. (3.30) 

This is also shown in the right-panel of fig. Differentiating ( p.26| ) with respect to time 
and repeating the steps in the previous section yields a simple expression for the rate of 
change of 99 in the minimum 



so the condition <C A'^^x^ implies that 

(1 + 4/3/7?) 



7 > 



6r/(l + 4/3r/) 



must be satisfied, which is true when x is affected by the DBI terms in the action. As with 
the case the of the offset potential, we find that the slow-roll parameter 

' 2(1 + 4/3/7?)^4y ^ ^2(1 + 4^/^) [2V J ^■^■■^"^^ 

is, as in the previous case, smaller than the standard DBI model when A > 1, leading to 
to a greater amount of efolds of inflation (see fig. Also, as in the example of the offset 
potential, one can see in fig. |^ a decrease in the boost factor of the DBI field. 
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4. Cosmological perturbations 



4.1 First order perturbation equations 



We now discuss the linear perturbations in our model. The model under consideration 
is part of a much larger class of multifield models with a Lagrangian of the form P = 
P{X^'^ ,(f>^), where X^"^ = —^d^(j)^d^(j)'^ is the kinetic term involving derivatives of the 
scalar fields (f)^ and cf)'^ . The equations of motion for the linear perturbations in the general 
case are calculated in [26|, but in order to render transparent the phenomenology that 
results from the presence of the coupling and to compare to the standard case, we shall 
perturb the field equations resulting from the action ( p. Ill) in the standard way. We 
decompose the fields tp and x i^ito a homogeneous and perturbed part 



V9(t,x) = ip{t) + (5v?(t,x), x(i,x) = x{t) +5x(i,x). 



We shall work with the Fourier components 5(p\^{i) and (5xk(i) so that spatial derivatives 
can be written V'^dip = —k'^Sip. Here, as below, we omit the subscript k to simplify the 
notation. In the longitudinal gauge and in the absence of anisotropic stress, the scalar 
perturbations of the FRW metric can be expressed as 



-(1 + 2^)dt^ + 0^(1 - 2'^)5,jdx'dx^ . 



(4.1) 



The perturbed Einstein equations give 



3H{^ + H^) + —-^ 



Sp, 



5q, 



Sp, 



(4.2) 
(4.3) 
(4.4) 



where 5p, Sq and 5p are perturbations of the total energy density, momentum potential 
and pressure respectively, given by 

Sq = -[ip{Sip)+jA''x{Sx)], (4.5) 
dp = [ipiSipY - ip^^] + A'-f^ixiSxY - x'*] + [U^ + 4^PDBi - ^^V'x'] ((5<^) + 



+ 



^4^^+ 1^4^(2 -7-^ -7) 



The perturbed field equations are 



k' 



(4.6) 
(4.7) 



{S^Y + 3H{d^Y + 

- (3A\3c;' - cj^)x{SxY + PA' 



+ [2U^ - 2PTDBI + pA\3c-^ - cf)i'] ^ - 49P* 
= 



(4.8) 



and 



(<^X)" + 



2,H + 2l5ip - 3- 



f: 



S 9 



s "XX 



-X-2X CsfxVx + 



f Cs 



f 



c-1 / JX 



fx 



f 



+ /3(3c2-l)x(M'+/3 



A^V^Csil + cl) + A2^c,(1 - csf + 2-^x 



- (3c^ + l)x* + 



A^V^Cs{l + cl) + A^^-^Cs{l-Csf + 2^X - 
(4.9) 



where we have rewritten factors of 7 using the expression for the sound speed for the DBI 
field = 7~^. In terms of the gauge invariant variables, 



Qx = (<5x) + I*, 



(4.10) 
(4.11) 
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the equations form a closed system, 



Qx + 



3H + 2/3^3 - 3— 

Cs 



(4.12) 



Qx + B^Q^ +{-ci + C^x Qx + CxvQv = (4.13) 



where we have defined the fohowing coeflacients 



B,r 



I3{3ci - 1) 
-A Cg Bip, 



X, 



(4.14) 
(4.15) 



a 



C 



XX 



a 



XV 



/5 ( I ) ^V'HSc. + 1)(1 - c-'f - ^^AT\IQ - 8ci' - 9cs + c;3) + 



2n a2^^X^ 



^'^^ fx „ /I -1\2/ -2 



+ 



H 



(4.16) 



rcs(l 



+ 



4^ P 

-lf3A^f-\3cs + l)il-c;'f 



c, - + 2cj' - 1) + 307^^2x0 - cj\l + c^)^^!!^ 



/ He, 



^-1^2 '^^X , "^^T.xVi 



2/72 



+ 



H 



+ 



^'X/X.. „^2 ^'X 



H 

X 

H p 



(4.17) 



+^A^1-Cs 



2H^ 
H 



f He 
fx\ 

X 

A^crVi + c 



-X - IcsPA~^x'Vt,x + 



+(l + c,)/ 



-1 I fx 

f 



,x 



+ lA\'c-Hl + e'g)- 



2„-l.. , „2^xV^X^(^^^2)^^3^-2^^_^^^ 



2/3 



4/^2 



H 



(4.18) 



2/3^r^(l-c.)2 



2H^ 



+ 



3V3X 



2if2 



+ 2/3csyl-VT,x + Cs- 



H 



+ 



H 



(4.19) 



In these expressions, Vt(¥',x) = U{ip) + A'^{ip)V{x) is the total potential. Comparing 



the form of these equations to |22|, where the authors treated the case of two coupled 
canonical fields, there are marked similarities. The principle diff'erences are due to factors 
of Cs that arise due to the non-canonical kinetic term of the DBI field. Let us define the 
auxiliary fields and u^, by 



'"'xQxj 



"X ~ ' x^x 



(4.20) 



with = aAc. 



'3/2 



and Tip = a, which correspond to the perturbation variables used in the 
analysis of the perturbations in single-field canonical and DBI infiation respectively. The 
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Figure 8: Perturbation quantities for tlie exponential case with 7^1 shown in fig. |^. The upper 
panels show the evolution of the amplitude of the variables (solid) and (dashed) for the 
first run (top-left) with Q"" = and the second run (top-right) with Q"' — for the pivot scale 
with /3 = 1. (Perturbations are normalised so that |(5™'| = 0.1 for the first run.) The centre plots 
show the evolution of the comoving curvature perturbation (centre-left) and the non-adiabatic 
contribution to its time derivative Hj (2if)(5pnad for the first (blue, dashed) and second (red, solid) 
runs in this case. Vertical dotted lines indicate the sound and curvature horizons for this mode. 
The lower panels show the variation in the resulting value of the spectral index rig (bottom-left) and 
power spectrum amplitude Pamp (bottom-right) with the coupling /3. The shaded regions shows 
observational values of these quantities (WMAP-I- BAO-l-HO) at 68% c.l il. 



equations of motion ( 4.12| ) and (4.12) can then be written 



^^x + 



T 

^x 
r' 



0, 



0, 



(4.21) 
(4.22) 



with B = r^B^. In this form the symmetry between the equations is more clear: each field 
perturbation exchanges energy with the other via the friction term B as well as the mass 
terms. The obvious complexity of these equations makes an analytical treatment difficult, 
so in the following section we shall integrate the perturbation equations numerically. 

4.2 The evolution of the perturbations 

In order to integrate the perturbation equations, it is important to specify the initial 
conditions. In the case of single-field k-inflation, the approach introduced by Garriga and 
Mukhanov |11] is to construct the auxiliary variable v = zlZ , with z = a{p + p)^^"^ /(cgH). 
One then quantises u, substituting u oc z/^e*^" to obtain, 



+ 



elk' 



z dr^ 



Vk = 0. 



(4.23) 



-20- 



Assuming the quantum fluctuations are in the trivial vacuum state, the solutions for fe^c^ ^ 
a^H^ are oscillatory as in the standard case, but with a frequency reduced by a factor of 
Cs that also affects the normalisation: 

i^k = -^=e-'^''^\ (4.24) 

The above procedure can be applied in the single field case, but when there are two 
fields present, one has to be to careful when choosing initial conditions for the numerical 
run. As in the standard two-field case we can neglect correlations between the perturbations 
for modes well within the horizon, so the short wavelength solutions for the auxiliary 
variables v^p and v-^ are 

so the relative normalisation of the two perturbations is dependent on the sound speed of 
the DBI field. However, in order to avoid introducing correlations between the fields at 
early times when integrating numerically (in which case the Bunch-Davies vacua would be 
inappropriate) the perturbations must be integrated twice, and the results combined to 
form the curvature power spectrum (see Appendix for a more in-depth description of the 
numerical method). In the first run, shown in the upper-left panel of fig. ^, the scalar field 
perturbations are set in the Bunch-Davies vacuum with Q™^ = and in the second run 
(upper-right panel of fig. |8|) Q^^' = and the perturbations of the DBI field are set in the 
trivial vacuum. 

The comoving curvature perturbation is given by, 

n = -^- -^^6q (4.26) 

P + P 



where 5q is the momentum potential given in eq. ( |4.5| ) and p and p are the total background 
pressure and energy density respectively. In terms of the gauge invariant variables and 
this yields the simple expression, 

'^ = 7^A'PQ^ + ^'^xQx]- (4.27) 
{-2H) 

where {-2H) = (p^ + A^^x^. 

When the fields are minimally coupled and both fields are evolving, as in the standard 
DBI case, the sound speed Cg affects the evolution of the perturbations so that Qy^ freezes 
in at the sound horizon [kcg = aH) while Q<^ freezes in at the curvature horizon k = aH. 
However, in the non-minimally coupled case, discussed in Sec. |3]^, the homogeneous field 
if is in the minimum of its effective potential and its perturbations are suppressed, so the 
primary contribution to the comoving curvature perturbation IZ (middle-left in fig. ^ 
is the second run, in which the Qy- is dominant. Thus in the coupled case, the curvature 
perturbation freezes in at the sound horizon of the DBI field, as in a single-field DBI system. 
The lower panels in fig. |8| show the variation in the resulting value of the spectral index 
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Figure 9: Background dynamics witli small 7 (exponential potential). In this case, as /3 increases 
from 0.5 to 4, the large coupling decreases the warp factor so x slow-rolls before reaching smaller 
values so 7 ss 1 for a considerable period. Quantities plotted are the same as fig. ^ Parameter 
values: A = 7.8 x 10", m = 2.8 x IQ-^, [/q = 7 x lO'^, n = 0.07, 77 = 2.8. 



Us (bottom-left) and power spectrum amplitude Pamp (bottom-right) with the coupling 
(3. In single field DBI, the power spectrum amplitude increases with the number of efolds 
of inflation. Here, the increase in (3 is accompanied by an increase in A^max, leading to a 
monotonic increase of Pamp with /3. The spectral index is much affected by the small sound 
speed of the DBI field in this example, which counteracts the red-tilt due to the deviation 
from a pure de Sitter stage, so that exhibits a slightly blue tilt. 

If the coupling is large, or the mass of the DBI field small, the value of 7 can be close 
to 1 when the modes that affect the CMB freeze in. An example of this is shown in fig. |9|, 
where increasing the value of /3 not only decreases 7 and increases A^'max but gives rise to 
a period of slow-roll at the start of the run. This can be seen by comparing the behaviour 
of the slow-roll parameter e and the boost factor 7 at the start and end of the run. At 
first 7 ~ 1; afterwards increasing to values < 20 when the effective DBI warp factor f 
becomes significant. In this case, the sound and curvature horizons are almost identical 
and the field perturbations behave as in the case of canonical fields, freezing in at A; = aH 
(see fig. 10). As in the previous case, the perturbations of the (p field are subdominant. 



and only give small corrections to TZ. It can be seen in fig. 10, that, due to this effect, the 
spectral index displays a red-tilt putting it within current observational limits. 

The rate of change of the curvature perturbation can be written. 



7^ 



H 



i-2H) 



P 

P 



(4.28) 
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Figure 10: The evolution of tlie perturbation quantities for tlie small 7 case plotted in fig. || for 
the pivot scale with /3 — 3.5. Quantities plotted are the same as in fig. ^. Note that, as 7 « f , the 
sound and curvature horizons (dotted, vertical lines) are indistinguishable and modes freeze-on as 
in slow-roll inflation, giving a red-tilted spectrum. 



where non-adiabatic pressure is defined by, 



P 



(4.29) 



The quantity 5pm = 5p — 3H6q is the gauge invariant comoving density perturbation. This 
appears on the RHS of the Poisson equation. 



o^Pmj 



(4.30) 



and thus decays in the long- wavelength/late-time limit. It is given explicitly by. 



6p„ 



+ 



2H 



+ [vt,^ - pA^rH^ - 57 + 7')] + 



[fx 



^T,x + |/3^'^(2-37 + 7') 



Qx- (4.31) 



In the standard two-field case, the non-adiabatic pressure perturbation consists of a 
contribution proportional to 5pm and another proportional to the entropy field perturbation 
6s. The former is neghgible at late times, but the latter can act as a source for the curvature 
perturbation even after the fluctuations grow larger than the curvature horizon. In order 
to obtain an explicit form for Jpnad) we can rewrite ( [4.31 ) as 



5p„ 



-^-§n + A\^'-j)xQ^ + 



^V^(7^-7)-f 



7^ 



+ [(3A^r\A-f - 37-^ - 7' - 2^3] + [i/^7'x' - 2(^%x)-] (4.32) 
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and combine with ( 4.2^ ) to get, 



P 



V-2iJ L J 



7 



— Qcr 



Here, we have introduced adiabatic and entropic fields, defined as 

1 



6s 



V-2H 



(4.33) 

(4.34) 
(4.35) 



in a similar manner to the coupled two-field system |21], so that the adiabatic field is 
related to the comoving curvature perturbation by 



7^ 



V-2H 

and the partial derivative of the total pressure with respect to the entropy field is, 
1 



(4.36) 



V-2H 
1 



-/3^5^-V/'(4-37-i-7))x]. 



1^347^/^(1 



7^/^(1 -7-')V 



(4.37) 



An interesting feature of eqn. ( [4.33| ) is that, unlike the canonical two-field case and the 
standard coupled case, this expression cannot be decomposed into purely entropic and 
adiabatic perturbation terms. The extra terms are expressed in the gauge-invariant com- 
bination Q^, 



-f^A-^fxQ^ + (73 - 7) 



1 '^X 



(4.38) 



The fact that this term appears like this suggests that the field decomposition is not 
appropriate in this case. We can also see this if we consider the form of the equation of 
motion for 



Qa 



1 



V-2H 

e 1 



+ 



«^ \/-2ij 



(4.39) 



where the ... represent terms in Qip and Qx as well as terms proportional to Q^p and Q-^ 
not proportional to 1? ja} . From this, we can see that the Qo- equation cannot be recast in 
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the form, x + c^^x + ... = 0, as a scalar perturbation can. The same is true of the entropy 
field. Even if the definition of the adiabatic field were changed to include a non-canonical 
kinetic term itself, i.e. 

crVP^Qa = vVP^Q^ + xVP^Qx^ (4.40) 

where Pzc^^ = —2H, no extra terms proportion to k'^/a^ would arise, so the original 
conclusion would stand. The description of non-adiabatic pressure perturbations in such 
a system as this could perhaps be better treated by the formalism introduced in p5|| , 
considering separate potential and kinetic fluids. However, as in this paper we focus on the 
coupled case in which the perturbations of ip are suppressed, non-adiabatic perturbations 
are negligible. We will thus discuss this issue in more detail in future work focusing on the 
non-minimally coupled case discussed briefly in Sec. 3.3 . 

In the non-minimally coupled case under consideration, we have found that the coupled 
DBI model is effectively a single-field system, and so the contribution of non-adiabatic fields 
is negligible. This can be seen in the middle-right panels of figs. ^ and in which the 
non-adiabatic pressure contribution to TZ, ^(5pnad); is shown to decrease exponentially 
as inflation proceeds. Thus, the curvature perturbation is constant to a high degree of 
accuracy on large scales. 



5. Discussion 



In this paper we have studied the consequences of coupling the DBI action to a canonical 
scalar field in a scalar-tensor theory. The motivation for this work is purely phenomenologi- 
cal: the DBI inflationary scenario has its motivation in string theory, but it is an interesting 
example of k-inflation in its own right. The additional scalar field might have its origin as a 
moduli field in a higher-dimensional setup, but we leave the explicit construction for future 
work (although, see on how DBI inflation could be related to scalar tensor theories). 

We have investigated the case in which the canonical scalar field is heavy during in- 
flation. We found this to be the case even for couplings of order a tenth of the strength of 
gravity (/3 ~ 0.1). Thus, we believe this to be a realistic scenario. With couplings larger 
than this, the additional scalar field is forced into the minimum of its effective potential. 
Thus, the parameters of the DBI field, which are dependent on the additional scalar field, 
vary during inflation. As a consequence, we find that the number of e-folds is extended and 
the boost factor is decreased (when compared to standard DBI inflation with the same bare 
parameters). The energy density during inflation comes mostly from the additional scalar 
field that sits in the minimum of the effective potential, which varies slowly as the DBI 
field evolves. But because the canonical field is heavy during inflation, its perturbations 
are negligible and only the DBI perturbations contribute to the curvature perturbation. 
Since the boost factor is smaller than in the standard case, the level of non-Gaussianity 
(which is proportional to 7^) is smaller than the standard DBI case. In fact, we found that 
the presence of the canonical scalar field alleviates the problems of the DBI infiationary 
scenario. 

We have considered two potentials for the canonical scalar field: a quadratic potential 
(with a minimum ^Pmin > 0) and an exponential potential. The difference between the 
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phenomenology of these two models is slight, mainly due to the presence of a minimum 
in the bare potential in the offset quadratic case, which affects the background dynamics 
towards the end of inflation, when the DBI contribution to the effective potential is less 
important. In the case of the exponential potential, the energy densities of the scalar fields 
(dominated by their potential energies) are proportional to each other (as in a tracker 
solution). In both cases we find that the behaviour of the perturbations is extremely 
similar, affected only indirectly by presence of the Lp field in its minimum and the details 
of the DBI dynamics; thus, we venture to suggest that our conclusions may hold for any 
choice of potential steep enough to allow ip to closely track the minimum of its effective 
potential. 

Finally, we have investigated two possible couplings. In the minimally coupled case we 
find find that if one fine-tunes the mass scales so that both fields contribute roughly equally 
to the total energy density for a considerable number of efolds, the dynamics of the DBI field 
can be significantly affected by the additional contribution to the Hubble damping, leading 
to smaller values of 7 (relative to the DBI dominant case). As both fields are evolving, the 
interplay between the field perturbations suggests that isocurvature fluctuations could be 
important. However, the minimally coupled case is extremely sensitive to small changes 
in the field masses and the initial conditions of both fields. In the non-minimally coupled 
case, we find that the coupled DBI model is effectively a single-field system, and so the 
contribution of non-adiabatic fields is negligible. Moreover, the curvature perturbation is 
constant to a high degree of accuracy on large scales. 

We have not discussed in detail the range of allowed parameter for the model but will 
study this question in future work. A range of parameters affect the prediction for the 
spectral index (for both the exponential and the quadratic potential) and we are currently 
scanning the parameter space numerically. Much has to be learned in this model, such as 
the decoupling limit /3 — )• 0. In this case, the canonical scalar field will become light and 
isocurvature modes become important. We expect that the signatures of non-gaussianity 
will be very different. 
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A. Numerical method 

In this short appendix, we devote a few words to an explanation of the general numeri- 
cal method used to integrate the perturbations. The pivot scale for the CMB (fcpivot = 
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0.002Mpc~^ at which the power spectrum ampHtude is calculated is identified with 
the mode that crosses the curvature horizon (k = aH) 55 efolds before the end of infla- 
tion (cf. |5^), when the slow-roll parameter e ceases to be smaller than 1. The spectral 
index is calculated by integrating the perturbation equations for /c-values in the range 
k G [/cpivot X 10^^, /cpivot X 10^], starting when kcg = bOaH (so the wavelength of the per- 
turbation is much smaller than the sound horizon of the DBI field) and ending 45 efolds 
from the end of inflation, at which time the modes are beyond the horizon and have frozen 
in. We do not integrate the equations ( |4.12 ) and (f4.13| ) with respect to conformal time 



or = logo, as is commonly done, but instead use logarithmic time x = logt, which is 
efficient in dealing with both the background and perturbation variables. In this case the 
oscillatory initial conditions ( [4. 251) for and are 

where L is a normalisation factor chosen so that Q^" = 0.1. Representing d/dx by ' and 
writing % = a' /a the derivatives are 

g; = -{n + Q^, Q'^ = -[n + p^' + i^cst^ Q^, (A.2) 

where as in the previous expressions, background quantities are evaluated at tini = exp(3;ini). 

As in the standard two-field case (cf .[^) in order to make sure that the modes 
are uncorrelated deep within the horizon, we perform two numerical runs, the In the 
first run is set in the Bunch-Davies vacuum ( |A.1| ) with Q™' = and in the second 
run Q™' = and the perturbations of the DBI field are set in the (quasi) Bunch-Davies 
vacuum. Evaluating the curvature perturbation ( 4.27] ) for each run yields the solutions TZi 



and 7^2) which can be combined to form the curvature power spectrum 

VTz = ^i\ni\^ + \n2\'). (A.3) 

The power spectrum amplitude (for the pivot mode fcpivot) and the spectral index can then 
be calculated at the end of the numerical run. 
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